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Abstract 

We compute the internal energy of different Ising type models, both long-range and 
short-range, under Tsallis statistics using the microcanonical and the canonical 
ensembles and we discuss under which conditions both ensembles give equivalent 
results. 
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1 Introduction 

In this paper we study different Ising type systems in the microcanonical and 
the canonical ensembles under the generalized statistics introduced by Tsallis 
[1]. This is based upon the following alternative expression for the entropy: 

i - y ■ v q 

s q = — (1) 



q-L 

depending on a parameter q. A probability Pi is assigned to each of the 
i = 1, . . . , W microscopic configurations. The set of probabilities {p{\ must 
be obtained maximizing S q under the appropriate constraints. The physical 
observables associated to a microscopic dynamical function O are computed 
as statistical averages (O) obtained using the rule [2]: (O) = J^iOiPi with 
Pi = vi/HjVj an d Oi is the value of O at the configuration i whose proba- 
bility is pi. Different statistical ensembles give rise to different probabilities pi 
and we address in this paper the question of equivalence between the micro- 
canonical (fixed energy) and the canonical (fixed average energy) ensembles. 
To be more specific we consider Ising type systems with Hamiltonian: 
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Where s n — ±1, and the indexes n,m = 1, . . . , N = L d run over the distinct 
pairs of sites of a regular <i- dimensional lattice of lattice spacing equal to 1 and 
with periodic boundary conditions; r nm is the distance between the spins n and 
m, and the parameter a sets the interaction range. The energy of configuration 
i is Ei (notice that the ground state has zero energy). The usual short range 
Ising model is recovered in the limit a — > oo where l/r° m — > unless n and 
m are nearest neighbors with r nm = 1. For ct < d the average energy per 
particle diverges in the thermodynamic limit and the system is said to be 
non-extensive. More precisely, a convenient scale for the average energy per 
spin in a finite system of size iV is given by [3]: E/N ~ TV = — 1 _ a J^ d ■ We 
see that for a > d the average energy per spin scales as a constant in the limit 
of large N, whereas for a < d, it grows with the system size, a non-extensive 
behavior. Throughout this paper we will be considering the following cases: 

(a) a — oo, d — 1, the short range one dimensional Ising model. 

(b) a — oo, d — 2, the short range two dimensional Ising model. 

(c) a = 0.8, d — 1, a non extensive, long range, one dimensional Ising model. 

(d) a = 0, the non extensive infinite range Ising model, whose properties are 
basically independent of the spatial dimension d. 

For each of these cases we will compute the internal energy E(T, N) (the aver- 
age value of the Hamiltonian (H)) as a function of the temperature T and the 
number of particles N, in the microcanonical and the canonical ensembles. We 
will compare the results obtained in both ensembles using the standard defi- 
nition of temperature as well as a recent proposal for a physical temperature 
[4]- 



2 The microcanonical ensemble 



The microcanonical ensemble is defined by fixing the energy E and setting 
Pi = to those configurations whose energy is not equal to E. The maxi- 
mization problem for the original entropic form S q given by Eq.(l) with the 
constraint of given energy E, and the normalization condition J2iPi = 1, has 
the solution of equiprobability: 

'q(e,n)-\ e l = E, 

Pi = \ (3) 
0, otherwise. 

Q(E, N) is the number of configurations with energy E for a system with N 
particles. Using this solution, the entropy as a function of the energy is given 
byS q (E,N) = (Q(E, N) 1 ~ g - !)/(! - q). Finally, the temperature is defined 
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by the thermodynamic relation ± = || = Q(E, N)~ q m( ^ N \ This relation 
allows us to plot the energy as a function of the temperature, E(T, N). 

We see that an essential ingredient in the microcanonical ensemble is Q(E, N), 
the number of configurations having energy E for a system with iV particles. 
The exact form of this function depends on the interaction range parameter a 
and on the spatial dimension d, although the total number of configurations 
is Y^e Q(E, N) = 2 N independent of a and d. We now specify this function in 
each of the four cases of interest: 

(a) a — oo, d — 1. The possible energy levels are E k — 4k, k — 0, . . . , N/2 
(assuming that N is an even number) and their degeneracy is easily computed 
as: n(E k ,N) = 2(Q. 

(b) a = oo, d = 2. The possible energy values are E k = 4k, k = 0, . . . , N. 
No closed exact expression is known for the function Q(E,N). However, the 
exact solution for the partition function of the L x M Ising model [5], has 
allowed Beale [6] to write a Mathematica program that can actually compute 
Q(E, N) for moderate values of N. Using this program we have generated the 
exact values of Q(E, N) for N = 32 2 or smaller. 

(c) a = 0.8, d — 1. This is the more complicated case, because no analytical 
expression or exact numerical values is available. For small values of N, up 
to iV = 34 we have made a complete enumeration of the W = 2 N configu- 
rations and their energies £j. For larger sizes, up to iV = 3000, the values of 
Q(E,N) have been obtained by using a numerical sampling method known as 
Histogram by Overlapping Windows (HOW) [7]. Details of the implementation 
of the method for this particular problem can be found in [8]. 

(d) a = 0. This case is equivalent to consider the Bragg-Williams approx- 
imation to the solution of the Ising model [9]. The energy levels are E k = 
2k(N — k), k = 0, . . . , N/2 (assuming again that iV is an even number) and 
the number of states is: Q(E k ,N) = 2(f) for k = 0,1,..., N/2 - 1 and 

n(E k ,N) = (^ /2 )foik = N/2. 



3 The canonical ensemble 



In the canonical ensemble we fix a value for the average energy E = (Ti). Max- 
imization of the entropy, Eq.(l), under this constraint and the normalization 
condition J2iPi = 1, leads to the following solution for the probabilities Pi [2]: 



Pr = I 



[l-(l-g) e< /T']^ 
£.[1-(1- ? ) £ ,7T'] 

o, 



J— 1 



1 - (1 - q)e l /T' > 0, 
otherwise. 



(4) 



Where we have defined 
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(5) 



and T is the (inverse of the) Lagrange multiplier used to enforce the con- 
dition of fixed average energy in the maximization procedure. Finally, it is 
possible to show the validity of the the Legendre structure of the resulting 
thermodynamics formalism by proving the relation [10,2]: l/T — dS q /dE. 

Although Eqs.(4,5) form a closed set, it is very difficult to use them in that 
form because the number of terms in the sums is the number of microscopic 
configurations, or 2^, an extremely large number. It is possible to simplify 
the calculations by rewriting those equations in terms of the number of states 
Q(E,N). More precisely, one can replace the sum over microscopic configu- 
rations (with 2 N terms) by sum over energy levels (with an number of terms 
proportional to NN), ^ — ► Y,k^( E k,N). Using the function Q(E, N) it is 
then possible to compute numerically the different sums and hence to com- 
pute the internal energy (and other quantities of interest) within the canonical 
ensemble formalism [8]. 

4 Results 

In Fig.l we plot the internal energy versus the temperature, in the micro- 
canonical and the canonical ensembles obtained using the appropriate number 
of states Q(E,N) for each case (a-d). It is shown that, within the accuracy 
displayed in that figure and in all cases, the microcanonical and canonical 
results agree for q < 1 but disagree for q > 1. 

Recently [4] it has been suggested that an appropriate definition for a physical 
temperature is 

^(i+v s «)(ip <«> 

In Fig. 2 we plot the internal energy as a function of this physical temperature 
in the four cases of interest (a-d). We can see from this figure that the results 
of the microcanonical ensemble collapse in a single curve for all values of 
q, i.e. the dependence on q disappears when using the physical temperature 
in the microcanonical ensemble. This collapse is in agreement with a recent 
theoretical analysis of the definition of physical temperature [11]. The results 
of the canonical ensemble show some zig-zags which should be repaired by 
using a Maxwell-type construction (although we have not been able to find 
an unambiguous way of performing such a construction). Those zig-zags have 
their origin in the necessary transformation to go from the parameter T' to 
the temperature T as given by Eq. (5), and a plot of T versus T' already shows 
that behavior. The mathematical origin of the zig-zags lies in the condition 
1 — (1 — q)ei/T' > of Eq. 4: when we sweep along the values of T", the 
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Fig. 1. Internal energy £ as a function of the temperature T for the Ising model 
defined by (2) in the cases (a) a = oo, d = 1, (b) a = oo, d = 2, (c) a = 0.8, d = 1, 
(d) a = 0. The number of spins is iV = 1000. Solid lines are the results obtained in 
the microcanonical ensemble, while symbols correspond to the canonical ensemble. 
From left to right curves, the values of q are q = 0.5, 0.9, 1, 1.1, 1.5. 



number of different energy levels with a non-zero probability increases in a 
discrete way, thus producing steps or zig-zags. The number of these steps is 
proportional to iViV, and in the thermodynamic limit we will get an almost 
continuous set of energy levels. These zig-zags are more pronounced the larger 
the interaction range. Except for this zig-zag behavior, we see again that 
the microcanonical and the canonical ensemble results agree for q < 1 and 
disagree for q > 1. The same conclusions are reached when studying the 
system magnetization instead of the internal energy. A more detailed study, 
beyond the scope of this paper, should consider also other properties such the 
heat capacity or the magnetic susceptibility. 

Notice that the ultimate reason for not having equivalence between the two 
ensembles for q > 1, is that fluctuations of the energy in the canonical en- 
semble can not be neglected. We have checked that this is indeed the case 
by computing the energy fluctuations cr{H) = \J (H 2 ) — (H) 2 as a function of 
the system size. We have checked that the fluctuations, when normalized by 
the scale of energy, NN, do not decay to zero for increasing N in the range 
of temperatures for which the microcanonical and canonical ensemble do not 
agree. For q < 1 fluctuations do decay to zero with the system size in all the 
temperature range. 
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Fig. 2. Internal energy £ as a function of the physical temperature T p h ys . Same 
cases and symbol meanings as in Fig.l. 
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